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ABSTRACT 

Aims. If a significant fraction of binary star systems spent some time as an inclined triple, either during their formation process or as 
the outcome of several close dynamical encounters in a crowded stellar environemnt, then the number of planets in binaries would be 
signifcantly lower than that around single stars. The stellar chaotic phase preceding the instability of the triple and the wide oscillations 
(-H ' in eccentricity and inclination of the companion star due the high mutual inclination between the companion and the singleton would 

quickly eject planets orbiting the binary in S~type orbits. 
I ■ Mefhods. We perform numerical simulations of the dynamical evolution of hierarchical triple star systems with planets hosted around 

O ' the primary star of the inner binary. Different values of mutual inclination, binary separation and singleton initial semimajor axis are 

explored in a statistical way. 

C/5 , Results. We find that a significant mutual inclination i„, between the singleton and the binary is a key factor for instability of the 

. planetary system. When /,„ is larger than ~ 40° the fraction of planets in the binary surviving the chaotic phase of the triple declines 

dramatically. The combination of eccentricity and inclination oscillations of the binary companion induced by the secular perturbations 
of the singleton and the sequence of close encounters preceding the ejection of one star fully destabilize a planetary system extending 
, beyond 1 AU from the star. For !,„ around 90° the percentage of surviving planets is lower than 20% for all binaries with a semimajor 

J!L . axis smaller than 200 AU. 

iir ' Conclusions. The frequency of planets in binaries with low separation may be strongly reduced by the residence of the pair in the 

' past in a temporary inclined hierarchical triple. 

' Key words. Planets and satellites: formation; Methods: N-body simulations; Celestial mechanics 

o 

C 1. Introduction bits around the baricenter of a binary system. A hierarchical 

triple can become unstable after some time, depending on its ini- 
L: Over 65 percent of the main sequence stars in the solar tial orbital parameters, leading to the disintegration of the triple 
. ^ neighborhood are members of binary or multiple star systems feggleton and KiselevaLfToollKiseleva et allfTOQfih . This disin- 
(Duquennoy and Mayor, 1991). As a consequence, answering tegration occurs via a phase of chaotic evolution whose outcome 
5-^ ; the question of whether planets can form and persist near one of t^e ejection of one of the three stars (typically the least massive 
. the stars in a binary has far-reaching implications for the overall bo^y) on an unbound trajectory. The other two stars, members of 
frequency of planetary systems. Studies on the long term sta- j^e original binary, are left in a more tightly bound binary. In a 
bihty of planets in binaries have shown that a planet cannot previous paper (Marzari and Barbieri, 2007, herein after MBl) 
be located too far away from the host star or its orbit will be showed that the orbital changes of the binary and the strong 
destabilized by the gravitational perturbations of the companion gravitational perturbations during the chaotic phase prior to the 
star. Holman and Wiegert (1999) found that the stable/unstable singleton ejection can influence the final configuration of a plan- 
boundary depends on the mass ratio and eccentricity of the bi- gtary system hosted by the primary star of the pair However, in 
nary, but for a wide range of parameters stable orbits may extend j^e context of near-coplanarity between the binary and the sin- 
well beyond one tenth of the binary semimajor axis. However, in ^^^^^^^ pi^nets can survive the triple decay in most cases and 
terms of the probability of finding a planet in binary systems the ^dapt to the new orbital parameters of the binary. The major ef- 
dynamical stabihty analysis is not exaustive since it does not take fg^t would be a significant change in the orbital configuration of 
into account the profund influence that stellar dynamic interac- t^e system after the triple instability with respect to the original 
tions may have had on flie early evolution of a planetary system configuration, as an outcome of the planet formation process, 
in a binary. ^^^^ paper we consider the dynamic effects of the decay 

It has been suggested that most binaries originate from the of inclined hierarchical triples on planetary systems. In partic- 
decay of multiple systems (Reipurth, 2000; Larson 1972, 1995, uiar, we will focus on planet survival during the unstable triple 
2001 ; Kroupa 1995). The most common configuration among configuration. At present, determinations of the mutual inclina- 
multiple systems is the hierarchical triple, where a singleton or- tjong two orbits in hierarchical triple stellar systems are 

available onl y for a very limited number of cases and are often 
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of 20 triple star systems have an inclination e xceeding 15° and 
are not coplanar.' Sterzik and TokovininI ( |2002|) analysed a differ- 
ent set of 22 visual triples finding an average mutual inclination 
of 79°. However, in both the studies the mutual inclination was 
derived from incomplete observational data. To compute unam- 
biguous mutual inclinations for triples, both radial velocities and 
visual orbits are required for the inner and outer system. So far, 
only six nearby systems have been observed with both methods 
and h ave direct and precise measured orbits (iMuterspaugh et al 
l2006t) . The values of mutual inclination for these systems range 
from 24° to 132° but the sample is too small to give hints on 
the real distribution of inclinations among triples. The mutual 
inclination of triples may either be primordial and related to the 
formation process of the triple by fragmentation of a molecular 
cloud or it may form at later times because of dynamical inter- 
actions, like encounters, between single stars and binaries in a 
dense cluster-like environment. In the latter case the inclination 
is due to the encounter geometry between the binary and the sin- 
gle star and should be randomly distributed. Any deviation from 
randomly oriented orbits may be an important indication of the 
relative importance of the two formation mechanisms. Assuming 
that planets can form in the binary before the bound hierarchical 
triple becomes unstable, the dynamical interactions between the 
stars during the chaotic phase can strongly affect the stability of 
the planetary system. 

We can envision two different scenarios for planet forma- 
tion and subsequent de stabilization within an inclined hierarchi- 
cal triple: 

- A primordial binary star system forms in a star cluster and 
planets accumulate from a circum stellar disk a r ound the 
main star either b y core-accre tion dPollack et all 1 19961) or 
by disk instability (jBosi,[l997'). The existence of a few gas 
giant planets in binary systems with separation of a few 
tens of AU (y Cephei and GL 86) suggests that the per- 
turbations of the companion star are not st rong enough to 
prevent the f ormat ion of planets in binaries dThebault et al.L 
l2006t iBossl 12006"). Successively, a temporary hierarchical 
triple builds up because of dynamical interactions between 
the primordial binary, with pla nets, and a passing by single- 
ton or binary dFord et all l2000i) . In a dense stellar environ- 
ment with a large abundanc e of binaries this is a frequent 
event ('McMill an et ai[T99Ih . In the presence of a significant 
inclination between the singleton and the binary orbit, the 
planetary system in the binary is strongly destabilized during 
the transitional triple state by both the secular perturbations 
of the singleton and the frequent stellar encounters during 
the chaotic phase preceding the break up of the triple into a 
binary and singleton. The destruction of the triple may occur 
either because it is unstable and after a short timescale the 
singleton escapes or because of an encounter with other ob- 
jects, single stars or binaries. After the ejection of the single- 
ton in a hyperbolic orbit, the primordial binary has different 
orbital parameters but it is also depleted, in most cases, of its 
original planetary system. 

- A primordial inclined hierarchical triple forms by fragmen- 
tation of a single, rotating, dense molecular cloud (Boss, 
[T98§). Planets can grow on S-type orbits around the primar y 
star of the binary, possibly by disk instability (iBossl[T997l) . 
Self-gravitating density clumps can contract into planets in 
only a few hundred years. It is unrealistic to expect that plan- 
ets can form by core-accretion in this scenario because of the 
strong secular perturbations that the singleton would apply 
on a long timescale on the secondary star and then, indirectly. 



to a putative planetesimal disk around the binary (assuming 
that planetesimals could form in such a highly perturbed cir- 
cumstellar disk). Even in this case, the combined destabiliz- 
ing effects of secular perturbations and stellar encounters in 
the chaotic phase preceding the ejection of the singleton in a 
hyperbolic orbit destabilize the planetary system around the 
primary. 

In this paper we numerically model the orbital evolution of 
planets in S-type orbits in a binary member of an inclined unsta- 
ble hierarchical triple. We consider different mutual inclinations 
between the binary and the singleton, while in MBl we simu- 
lated only the planar case. We find that, contrary to the low in- 
clined cases, only in a limited number of cases do planetary sys- 
tems extending beyond 1 AU survive after the chaotic phase of 
stellar encounters when the mutual inclination between singleton 
and companion is larger than ~ 40°. There are two mechanisms 
that, acting in synergy, destabilize planets around the primary 
star. 

- A large initial mutual inclination between the outer stars ex- 
cites consistent eccentricity-inclination oscillations of the 
binary companion with periods of the order of some thou- 
sand years and more. These oscill ations, well described by 
a quadrupole-le vel secular theory dMazeh & Shahanl 119791 : 
iFord et all 12000) . strongly affect the orbits of the planets 
around the primary, forcing most of them to leave the system 
on hyperbolic orbits. In the phase of high eccentricity, the 
companion moves clo ser to the planetary system, reducing 
the region of stability dHolman and Wiegerll 1 19991) and per- 
turbing the planetary orbits. After a few Kozai-cycles a large 
fraction of the planetary system is destroyed. For retrograde 
orbits the dynamics is more complex but the evolution is still 
characterized by wide oscillations of both eccentricity and 
inclination. This kind of Koza i mechanism in s t ellar tr iples is 
different to that described by iMalmberg et alj d2007[) where 
the secular interactions involved only the companion star of 
an isolated binary and the planets. 

- The second mechanism giving the 'coup de grace' to the 
planetary system is the sequence of close encounters be- 
tween the singleton and the companion star of the binary 
occurring during the chaotic phase preceding the disgrega- 
tion of the triple. Close encounters between stars on mutu- 
ally inclined orbits are much more effective in destabilizing 
the planets than in the low inclination case. 

Both the eccentricity oscillations of the secondary star and 
the close encounters conspire against the survival of a planetary 
system around the main star of the binary. If indeed binary stars 
are born as triple or higher multiplicity stellar systems or they are 
temporarily involved in unstable triples, their planetary systems 
would be fully destabilized in most of the inclined cases. The 
fraction of planets in binaries observed at present would than 
be lower than that around single stars. The number of binaries 
depleted of planets would depend on the distribution of mutual 
inclination in the primordial triples. If coplanarity dominates, 
as might be the case if the majority of temporary triples orig- 
inated direct ly from the fragmentation of an interstellar cloud 
dBossl Il988l) . then the binaries stripped of their planets would 
be a minority. On the other hand, if most of the triples formed 
by gravitational interactions in a dense stellar environment, the 
large mutual inclinations would lead to a strong planet depletion 
among the surviving binaries. In this scenario, it is an important 
observational challenge to increase the statistics of known triple 
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systems with unambigous determination of the mutual inclina- 
tion. At present only about 15% of planets have been found in 
multiple stellar systems. This is probably an observational selec- 
tion effect in favor of single target stars but it might also be a first 
indication that binary systems are depleted of planets by the past 
violent dynamical evolution of the stars. 

We will not explore in this paper the full complexity of the 
hierarchical triple dynamics as performd in Ford et al. (2000). 
We are interested on the consequences of the large variations of 
the star orbital elements on the planets and we perform statisti- 
cal numerical simulations giving the fraction of planets surviv- 
ing the chaotic phase of unstable triples. We also do not inves- 
tigate the planetary formation process in detail, but we assume 
that planets can form by either of the two mechanisms, core- 
accretion or gravitational instability. 

In Sec. 2. we describe the numerical model adopted for the 
numerical integration of the trajectories of the stars and planets. 
Sec. 3 is devoted to the statistical analysis of the survival of plan- 
ets in S-type orbits around the primary star. In Sec. 4 we present 
our conclusions. 

2. The model 

Our numerical model consists of 3 stars, two locked in a binary 
system and the third orbiting the barycenter of the pair A set of 
10 massless bodies started on circular orbits around the primary 
star simulate a planetary system that formed in the early phases 
of evolution of the binary. The semimajor axes of the test bodies 
are regularly spaced from 1 to 10 AU and the initial inclinations 
are all set to 0° with respect to the binary orbital plane. The tra- 
jectories of the stars and of the 'planets' are co mputed with the 
numerical integrator RADAU (lEverhartl Il985l) . It handles in a 
very precise manner close encounters between massive bodies 
and it does no t require a fix ed hierarchical structu re such as HJS 
jBeusl 12003) or SYMBA5 (iDuncan et al.LfT998h . 

To model the outcome of the triple instability in all possi- 
ble configurations is a difficult task since the parameter space to 
explore is wide. For this reason we select a limited number of 
parameters to be sampled while the others are left unchanged. 
To better compare our results with those presented in MBl, we 
adopt the same masses for the stars i.e. 1 and 0.4 solar masses 
for the binary, 0.4 solar mass for the singleton. An eccentricity 
of 0.2 is adopted for both the binary and the singleton, taking 
into account that the orbit of the singleton is defined with re- 
spect to the barycenter of the binary. The mutual inclination is 
sampled between 0° and 180° including in this way retrograde 
orbits of the singleton. For any value of the semimajor axis of 
the binary cii,, we sample different values of the semimajor axis 
of the singleton Oj and of the orbital angles other than those giv- 
ing the mutual inclination. For any set of {at,, as, /,„) we perform 
20 simulations with random initial orbital angles to increase the 
statistics on the star and planet dynamical behaviour. 

3. The dynamical sources of instability 

In this section we discuss in detail the two mechanisms leading 
to destabilization of a putative planetary system extending be- 
yond 1 AU around the primary star of a binary in an inclined 
temporary triple. In FiglT] we show the evolution of a model 
with ab = 70 AU, a., = 212 AU and initial mutual inclination 
!,„ = 90°. In this configuration, the critical semimajor axis for 
long-term stability of planetary orbits around the primary is, ac- 
cording to lHolman and WiegertI ( II999 ). around 21 AU. Our ini- 
tial planetary system, extending out to 10 AU, is well within the 



stability region. The perturbation of the singleton induces Kozai 
cycles on the binary companion that achieves an eccentricity of 
almost one over a timescale of 2.5 x lO** yrs. This behaviour is 
well described by quadrupole and octupole - level secular equa - 
tions described in Mazeh & ShahanI (Il979l) : iFord et al.1 (|2000|) . 
All the planets beyond 2 AU are ejected from the system after 
the first cycle, while that orbiting at 2 AU is destabilized after 
the second cycle. Starting from 1 x 10^ yrs the singleton and 
the companion star have mutual close encounters that quickly 
lead to the ejection of the last inner planet, that lived through the 
Kozai cycles of the companion. Finally, after about 5x10^ yrs, 
the outer star is ejected on a hyperbolic orbit and the the binary 
system is left with a smaller separation but no planets. 

This kind of behaviour, typical of systems with high mu- 
tual inclination /„,, places in jeopardy not only the stability of 
planets around the primary but also the possibility that they can 
form. According to iBossI (11997), several gaseous protoplanets 
can rapidly form by disk instability in a marginally gravitation- 
ally unstable protoplanetary disk. Within this scenario in a few 
hundreds years we might witness the formation of the unstable 
triple and of a planetary system made of gas giant planets around 
the primary before the Kozai cycle increases the eccentricity for 
the companion star On the other hand, core-accretion would not 
have enough time to accumulate a core by planetesimal accre- 
tion, and even planetesimals may have failed to form on such a 
short timescale. A protoplanetary disk around the primary star 
would be strongly perturbed and almost fully destroyed during 
the first Kozai cycle in eccentricity and inclination of the com- 
panion star. However, if the binary system was isolated during 
its formation and it became part of an unstable triple later on 
because of repeated stellar encounters in a dense star-forming 
region, then planets might have the time to grow even by core- 
accretion, before the onset of the strong perturbations related to 
the stellar interactions in the triple phase. 

Only systems with large values of /,„ are fully destabilized 
by the stellar perturbations of the triple. When the mutual incli- 
nation is lower than ~ 40°, close encounters between the stars, 
and the consequent impulsive changes of the orbital elements, 
are a source of instability for the planets but often not strong 
enough to destabilize the full planetary system. In Fig|2]we il- 
lustrate the evolution of a model with ah - 70 AU, Oj - 212 
AU, as in the previous case, but with a lower initial mutual in- 
clination /„, - 30°. The triple quickly becomes unstable and the 
singleton has frequent close approaches with the binary compan- 
ion marked by sudden steps in eccentricity and semimajor axis. 
The changes in the orbital elements of the companion leads to 
unstable planetary orbits as shown in Fig|2l However, contrary 
to the case shown in FiglT] the planetary system is not fully de- 
stroyed and planets within 5 AU of the star survive the chaotic 
phase. Further perturbations by the binary companion after the 
triple disruption do not destabilize the planetary survivors since 
they are well within the criti cal semimajor axis for stable orbits 
dHolman and WiegertilT999l) . 

If the companion and the singleton are on retrograde orbits, 
instabiUty builds up in a similar way. For mutual inclinations 
lower than ~ 140°, large amplitude oscillations of the eccentric- 
ity begin to destabilize the planetary system which is finally de- 
stroyed by the stellar encounters in the chaotic phase. Contrary 
to the prograde case, the oscillations of eccentricity and incli- 
nation are not in phase, as predicted by the quadrupole theory. 
Apparently, the two orbital parameters are no longer bound in 
an invariant and they evolve with independent frequencies. For 
inclinations in between ~ 140° and ~ 180° the oscillations in ec- 
centricity are moderate to low but some instability of planetary 
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Fig. 1. Orbital evolution of the planets around the primary star of 
the binary (plot a) under the perturbations of the companion star, 
in turn affected by the gravitational pull of the outer singleton 
star (plot b). The initial semimajor axis of the binary is 70 AU, 
the eccentricity of the binary 0.2, that of the singleton 0.2, and 
the mutual inclination /„, is set to 90°. 
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Fig. 3. Destabilization of planets around the primary when the 
singleton is on a retrograde orbit relative to the companion. The 
initial mutual inclination between the two outer stars is 150°. 
The AQ angle librates around 180° while the inclinations of the 
two planets have wide oscillations. Only two planets, the closer 
ones, survive after 1 Myr of evolution. 
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orbits is driven by the large inclination oscillations of the com- 
panion. In Figl3]we show the evolution of the inclinations of the 
singleton and companion star when the mutual initial inclination 
is ~ 150°. We plot the inclination of each individual star, referred 
to the initial plane of the binary, because this is the plane where 
the planets also begin to orbit the primary. The inclination of 
the companion becomes very high and becomes retrograde for a 
short while. The behaviour is characterized also by the libration 
of the angle Af22 - AQi (see Fig|3] lower panel) with the same 
frequency as the inclination oscillations. Most of the planetary 
orbits are destabilized during these large inclination excursions 
of the companion star and after 1 Myr only the two inner planets 
survive. The onset of the chaotic phase of the stars ejects finally 
also these two survivors. 
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Fig. 4. Fraction of binaries that, at the end of the temporary 
triple phase, retain at least one of the planets in a stable orbit 
around the primary vs. /„„ the initial mutual inclination between 
the singleton and the binary companion. 



4. Statistical outcome 

To test the chances of a planetary system in a binary to survive 
a period of stellar interactions typical of an unstable triple, we 
have run several models with the binary semimajor axis fixed to 
at, = 70 AU. The orbital eccentricities of the stars are both set 
to 0.2. In Fig|4]we plot the percentage Pj of dynamical systems 
that, at the end of the period as a hierarchical triple, retain at 
least one of the initial 10 planets vs. the initial mutual incli- 
nation between the two outer stars. This percentage is very high 
for low inclinations confirming the results presented in MB 1 for 
low-inclination systems, while it declines very quickly when the 
inclination approaches 90°. This is a consequence of both the 
Kozai cycle that pushes the binary companion closer to the plan- 
ets, and of the more complex orbital behaviour during close ap- 
proaches between the stars when their orbits are inclined. 

Retrograde orbits of the singleton also lead to fast instability 
of the planets when the mutual inclination is close to 90°. Wide 
oscillations of the eccentricity up to large values are observed, 
even if not related to the known Kozai type mechanism: there is 
no phasing between eccentricity and inclination. However, even 
in this case when the eccentricity is at its peak value most of the 
planets are destabilized. Only when the mutual inclination has 
values beyond ~ 140° the planets around the primary are partly 
spared by the oscillations in eccentricity of the binary compan- 
ion. However, as observed in the previous section, for mutual 
inclinations in the range ~ 140° - 180° large inclination oscilla- 
tions of the companion star destabilize planets even if to a lesser 
extent than the eccentricity oscillation. As a result, the percent- 
age of planets surviving the chaotic phase grows for inclinations 
larger than 90° but it does not return to 100%, halting at about 
30%. 

If we increase the semimajor axis of the binary ai,, the frac- 
tion of systems with surviving planets increases in an almost lin- 
ear way. In Fig|5]we show the fraction of systems retaining plan- 
ets vs. Ub for the worst case, i.e. with mutual inclination equal to 
90°. The triple instability is a mechanism that easily destroys 
planetary systems of close binaries while it is less effective for 
wide binaries. For larger values of ai, the planetary systems that 
survive are also more extended in semimajor axis. In most cases 



Singleton semimajor axis (AU) 

Fig. 5. Fraction Ps of binaries that, at the end of the temporary 
triple phase, retain at least one of the planets in a stable orbit 
around the primary vs. at,, the binary semimajor axis. The mutual 
inclination /„, is set to 90°, the worst case for planet stability. 

for at, = 250 AU all the planets up to ap - 10 AU survive the 
stellai- chaotic phase. 

5. Conclusions 

The fraction of binary systems hosting planets in S-type orbits 
can be lower than expected. If the binary is part of a crowded 
stellar environment, encounters with other stellar objects can 
lead to the formation of a transitional triple with large mutual in- 
clination between the singleton and the binary. The subsequent 
dynamic evolution of the triple, in particular the large oscilla- 
tions in eccentricity of the companion star in the binary and the 
chaotic evolution during the triple destruction, destabilize plan- 
etary orbits around the main star. Even if the binary was born as 
part of an unstable inclined triple, the planetary system is fated 
to be disrupted. 

Observing a binary system without planets in S-type orbits 
does not necessarily imply that the stars did not posses circum- 
stellar disks in their early phases or that planet formation did not 
occur The history of the binary and of its primordial environ- 
ment must be taken into account since in most cases it may be 
the cause of the absence of planets. Planet formation might be a 
very efficient process also in the presence of external perturba- 
tions, but the survival of planetary systems may be threatened by 
the binary dynamical history. 

Acknowledgements. We thank P. Eggleton for stimulating us to perform this 
work. 
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